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^ ^ Abstract. In this paper we present a new formula relating Stieltjes 

numbers 7^ and Laurent coefficinets i]^ of logarithmic derivative of the 

Qs ' Riemann's zeta function. Using it we derive an explicit formula for 

the oscillating part of Li's numbers A„ which are connected with the 
Riemann hypothesis. 

-)-> ■ 1. Introduction 

First recall some basic definitions and conventions. The Stieltjes constants 
7„ are essentially the coefficients in the Laurent expansion of the Riemann 
zeta function about s = 1: 



pg , zeT;a luncLion aooux s = i: 
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rj. (1.1) ^(, + i) = _ + ^^. 
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^O I (Sometimes another convention is adopted. It differs from the above by the 

factor (—1)" /n\) It is well known that zeta has a single simple pole at s = 1 



(•~^ . with residue 1 which is evident from (|l.lj) . 



Other useful coefficients t]^ are these which appear in the Laurent expan- 
sion of the logarithmic derivative of zeta about s = 1: 



(1-2) -^-{s + l) = -+^VnS'' 

K> ' It can be shown that IT] 



(1.3) ,„;=uH:,i^ vfl^i^-fi^:!^:' 
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where A{k) denotes the von Mangoldt function (see e.g. |^, p. 50). It is 
related to the Riemann zeta function by 

C^ (.) _ A A (n) 
C(.) ^^ ns 

where 3f?s > 1. It may be shown that A (n) is zero unless n is a prime power: 
n = p^ , in which case A (n) is log(p) 

However, the expansions (|1.3() and (|1.4() are slowly convergent and are 
therefore useless in numerical computations. Effective numerical algorithms 
for calculating 7^ were given by Keiper ^ and Kreminski |Hj . An interesting 
recurrence relation involving both 7„ and r/^ has been recently discovered 
by Coffey |2j (written here in our convention (|1.H) and (|l.Hj) concerning 7^): 

n-l 

(1-5) 'n^ = -{n + l)-f^-^ri,,-f^_k_i 

k=0 

2. Derivation 
Here are several initial relations obtained directly using ()1.2[) : 
(2.1) 



yo72 - 473 

57o72 + 57i72 + 57o73 - ^74 



With growing n the number of terms in 7y„ increases as the so-called partition 
number. 

At the first glance there seem to be some regularities in (|2.1j) . Indeed, 
careful (and very tedious) inspection of the indices in 7s, their exponents, 
as well as the numerical coefficients reveals several simple rules: 

(1) Each 7y„ is a linear combination of products of powers of Stieltjes 
numbers with some numerical coefficient: ibA7j'7''...7^. 

(2) Every such term appears only once. More precisely, ±^7^*7 .^...7^ 
appears in ry„ with n = {i + l)ki + {j + l)kj + ...{m + l)km- 

(3) The sign of each term is determined by the sum of all exponents kf. 

n 

it is equal to (—1)^, where p = Y2 ^i- 

i=0 

(4) The numerical coefficients A which appear in rj^ may be obtained as 

follows. Let us discard for a while in each term power of 79, i.e. 7g'', 
if there is one. It turns out that A contains two factors. The first is 

n 



Vo = 


-7o 


Vi = 


+7o - 271 


V2 = 


-7o + 37o7i - 372 


Vs = 


+7o - 47o7i + 27? + 


rii = 


-7o + 57o7i - 57o7? 



and the other is 
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q-2 

n {n-i) 

i=s+l 



where s is the sum of all indices of 7 in a given product and q is the 
number of all 7s in the same product (in both cases 7s are counted 
according to their multiplicity and in both cases we neglect 7q), i.e. 

s = iki+ jkj + ... + mkm 

Q ^ Ki ~r Kj ~r ■■■ ~r K-m 

The above product may be further simplified: 



r (n — s) 



Jii^^ '>- ^ ') ris-n+l)-r{n-s-q + l) 

where the last equality stems from the properties of the Pochham- 
mer's symbol. According to (2) we have n — s — q + 1 = ko + 1, 
and therefore T [n — s — q + 1) = k^l The last step is equivalent to 
restoring the powers of 7q in the general term ±747^*7-^. ..7^. 

Collecting all the above rules (1-4) together leads to the general expression 
for ?7 as a function of appropriate 7s: 



(2.2) 

where 
(2.3) 




p = p{ko,ki,...,kn) = y^fcj 

j=0 
n 

r = r{kQ,ki,...,kn) = '^{l + i)ki 



i=0 

and the sum (|2.2j) is performed over all combinations of integers ki = 0, 1, 2... 
satisfying the constraints (|2.3j) . This sum contains formally many terms, 
roughly 0{n'^), but, due to the Kronecker delta which picks out only appro- 
priate terms (i.e. when r = n) and cuts off all the others, the number of 
non-zero terms is much smaller. This neat and concise relation (|2.2j) is in 
fact pretty sophisticated. 
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Formula (|2.2j) may be inverted. Writing down several initial expressions 
(by solving the quasi- linear system (|2.1j) with respect to 7„) we get: 

(2y4) = -rio 

7i = 2 (+^0 ~ ^i) 

72 = ^ (-^0 + 3??o^i - 2r?2) 

73 = 77 (+^0 - 6r?o^i + Sr]l + 8r/o7?2 - 6??3) 



4! 
1 

5! 



74 = 77 (-^0 + 10??o??i - 15r/or?? - 20r?27?2 + 20r?i7?2 + SOryo^a " 24?74) 




Again, several regularities are evident. When collected together in the same 
way as before they give: 



(2.5) 



where r is given by the second relation (|2.3() . It should be emphasized 
that formulas 1)2. 5|) and 1)2. 2|) . as well as 1)3. 2|) below, may be effectively 
implemented e.g. using Mathematica symbolic package (it is not trivial 
since the number of summations and products is variable): 

Table[{Subscript [\ [Eta] ,n-l] , 

Factor[n*Sum[r[Sum[Subscript[k,i],{i,0,n}]]* 

KroneckerDelta[n,Sum[(l+i)*Subscript[k,i],{i,0,n}]]* 

Product [(-Subscript [\ [Gamma] ,i])" Subscript [k,i] /Subscript [k,i]!,{i,0,n}], 

Evaluate[Apply [Sequence, Table[{Subscript[k,j],0,n},{j,0,n}]]] 

]]},{n,l,5}] 

where the Euler gamma function T above is slightly modified: 

r[n_] := Gamma[If[n == 0, 1, n]] 

3. Applications 

It was shown elsewhere by the author [7\ that the behavior of certain 

numbers A„ with growing n is crucial for the Riemann hypothesis to be 
true. Now using 

(3-1) ^" = -E(-)'^.-i 
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together with (|2.2|) we get 



(3.2) 




which exphcitly expresses the oscillating part of Li's numbers using Stieltjes 
constants. Similarly as in (|2.2|) . due to the binomial coefficient there are in 
fact many redundant terms in 1)3. 2|) which are identically zero, i.e. always 
when r is greater than n. Those which are non-zero (r < n) give: 



(3.3) Ai 


= 7o 


A2 


= 270-70 + 271 


As 


= 37o - 37^ + 7g + 671 - 37o7i + 372 


A4 


= 470 - 67^ + 47|i - 7o + 1271 - 127o7i + 47o7i 




+1272 - 47o72 + 473 



2ii 



With growing n the number of terms in A„ increases as the summatory 
function for partition number, i.e. faster than in the case of 1)2. 2|1 . 



4. Discussion 

The distribution of values for various terms which contribute to (|3.3|) is 
very nonuniform. When n ^ 00 the majority of values is concentrated near 
zero. In Figure 1 there are eight histograms for n = 3, 4, ..., 10. The shape of 
these histograms is roughly symmetrical about zero but it is just the slight 

departure from symmetry which contributes to A„. 
(Figure 1 about here) 
The fundamental problem, equivalent to the Riemann hypothesis itself is 

whether the oscillating part of Li's numbers A„ is bounded by the trend A„ 
which grows asymptotically as: 



A. 



1 



(1 + nlnn) + en 



where (cf. [HI, 0): 



-(7 

2 ^' 



1 - In 27r) 
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Figure 1. 

Distribution of terms in formula (3.2) for 72=3,4,..., 10. These terms sum up to 
the oscillating part of Xn. As usual, the base of each rectangle is equal to the 
width of the class interval and its height is proportional to the number of data 
values in the class. The base is always equal to 1/5. Note changes in the hori- 
zontal scale when n increases. 



